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ABSTRACT 

We explore the effects of stellar population models on estimating star formation histories, ages and 
masses of high redshift galaxies. The focus is on the Thermally-Pulsing Asymptotic Giant Branch 
(TP-AGB) phase of stellar evolution, whose treatment is a source of major discrepancy among different 
evolutionary population synthesis. In particular, besides the models usually adopted in the literature, 
we use models (by Maraston 2005), in which the contribution of the TP-AGB phase is calibrated 
with local stellar populations and is the dominant source of bolometric and near-IR energy for stellar 
populations in the age range 0.2 to 2 Gyr. These models also have an underlying different treatment 
of convective overshooting and Red Giant Branch stars. For our experiment we use a sample of 
high-z (1.4^2^2.5) galaxies in the Hubble Ultra Deep Field held to be mostly in passive evolution, 
with low-resolution UV-spectroscopy and spectroscopic redshifts from GRAPES, and Spitzer IRAC 
and MIPS photometry from the Great Observatories Origins Deep Survey. We choose these galaxies 
because their mid-UV spectra exhibit features typical of A- or F-type stars, therefore TP-AGB stars 
ought to be expected in post-Main Sequence. We find that indeed the TP-AGB phase plays a key 
role in the interpretation of Spitzer data for high-z galaxies, when the rest-frame near-IR is sampled. 
When fitting without dust reddening, the models with the empirically-calibrated TP-AGB phase 
always reproduce better the observed spectral energy distributions (SEDs), in terms of a considerably 
smaller x^ ■ Allowing for dust reddening improves the fits with literature models in some cases. In 
both cases, the results from Maraston models imply younger ages by factors up to 6 and lower stellar 
masses (by ~ 60% on average) . The observed strengths of the Mgyy spectral feature compare better 
to the predicted ones in the case of the Maraston models, implying a better overall consistency of SED 
fitting. Finally, we find that photometric redshifts improve significantly using these models on the 
SEDs extending over the IRAC bands. These results are primarily the consequence of the treatment of 
the TP-AGB phase in the Maraston models, which produces models with redder rest-frame optical to 
near-IR colors. This work provides the first direct evidence of TP-AGB stars in the primeval Universe. 

Subject headings: stars: AGB and post-AGB — galaxies: evolution — galaxies: formation — galaxies: 
high-redshift 



1. INTRODUCTION 

Age-dating the stellar populations of galaxies provides 
us with a cosmic timescale which is independent of cos- 
mological models, and with a mean of reconstructing 
their star-formation history. Moreover, the same pop- 
ulation synthesis tools used for age dating provide also 
an estimate of the stellar mass of galaxies. Thus, deriv- 
ing star formation histories and masses for large num- 
bers of galaxies at various redshifts allow us to attempt 
an empirical reconstruction of the formation epochs and 
mass-assembly history of galaxies. Theoretical models 
for galaxy formation and evolution can then be adjusted 

^ Oxford University, Denys Wilkinson Building, Keble Road, 

OXl 3RH, Oxford, UK 

Electronic address: ma raston@astro.ox.ac.u k 

^ National Optical Astronomy Observatory, 950 N. Cherry Ave., 
Tucson, AZ, 85719 

^ INAF-Osservatorio Astronomico, Vicolo dell'Osservatorio 5, I- 
35122 Padova, Italy 

** INAF-Osservatorio Astrofisico di Arcetri, L.go E. Fermi 5, I- 
50125, Firenze, Italy 

^ Max-Planck- Institute for Astronomy, Konigstulil 17, 69117 
Heidelberg, Germany 

® Space Telescope Science Institute, 3700 San Martin Drive Bal- 
timore, MD 21218 USA 

'' Spitzer Fellow 



in order to comply with the galaxy evolution scenario 
emerging from the observations. 

Clearly, the evolutionary population synthesis tools 
play a critical role in this process, and we must pay as 
much attention as possible to ensure that they provide 
ages and masses as accurate as possible. In this paper we 
focus on the problem of age and mass determinations for 
galaxies which are dominated by stars in the age range 
0.2 ^ t/Gyr ^ 2 Gyr, because this is a range where differ- 
ent population synthesis tools give quite discrepant re- 
sults. Moreover, at redshifts beyond 2 (or 3) the universe 
itself is younger than ~ 3 (or ~ 2) Gyr, and galaxies at 
such high redshift are bound to contain stars younger 
than ^ 2 Gyr. Hence, the use of an accurate stellar pop- 
ulation tool is essential for the correct interpretation of 
the properties of high redshift galaxies, right at a cosmic 
time when different renditions of the hierarchical galaxy 
formation paradigm diverge dramatically (see e.g.. Fig. 
9 in Fonta na ct al. 2004). 

In stellar populations with ages 0.2^i/Gyr^2 Gyr, 
stars in the advanced evolutionary stage known as the 
thermally-pulsing asymptotic giant branch (TP-AGB) 
contribute up to r^ 40% of the bolometric light, and up to 
- 80% of the near-IR ( Maraston 199 8, hereafter M98). 
The first appearance of TP-AGB stars in the life of a 
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stella r population has been cal led the AGB-p/iase transi- 
tion IjR.enzini fc BuzzoiiiII198(tD . and its direct effect is a 
strong increase of the near-IR flux, as TP-AGB stars emit 
most of their radiation at A > Ifim. With the advent of 
the Spitzer Space Telescope (SST), the rest- frame near- 
IR is now accessible to very high-redshifts, which allows 
one an improved analysis of the spectral energy distribu- 
tion (SED) of high redshift galaxies, and makes it even 
more important than before to dispose of reliable stellar 
population tools. 

The following concrete example illustrates the case. 
For a sample of high-z galaxies with photometric red- 
shifts <z>~ 2.4 and with data extending from the rest- 
frame UV to near-IR (from SST) Yan et al. (2004) esti- 
mated an age of 1.5 to 3.5 Gyr for the dominant stellar 
population (~ 99% by mass), with the residual trace of 
young stars as due to a recent star-formation episode. 
The Yan et al. res ult was obtained u sing the synthetic 
population tool by iBruzual fc Chario t (2003, hereafter 
BC03). S ome of the same galaxies were re-analysed by 
iMarastonl (^2005, hereafter M05) using her own synthetic 
stellar populations, obtaining ages ^ 0.6 — 0.8 Gyr, right 
in the middle of the epoch when TP-AGB stars dominate. 
Moreover, the derived stellar masses are smaller when us- 
ing Maraston models ( <^ lO^^M© versus the typical value 
of ~ IO^Mq quoted by Yan et al. 2004). Whereas both 
results may be astrophysically plausible, worth empha- 
sising is that their implications concerning the formation 
redshift for the bulk of the stars and the galaxy mass 
assembly are vastly different. The difference is likely the 
result of the different weight given to the TP-AGB in the 
two sets of models, with the TP-AGB contribution being 
systematically lower in the BC03 models compared to the 
models of Maraston (2005). When using BC03 models a 
massive contribution of old stars is required in order to 
match the observed near-IR fluxes, and in the meantime 
a sprinkle of young stars is necessary to account for the 
rest-frame near-UV. On the contrary, a two-component 
stellar population is not required by the M05 models, as 
a single age (younger) population accounts at once for 
the SED from the near-UV to the near-IR. It is interest - 
ing to note in this context the early result of lLillvl l)1987f) . 
who found for a sample of z ~ 0.45 cluster galaxies that 
synthetic near-IR colours not including the AGB phase 
were not red enough to match the data at given optical 
colour. The problem was cured by adding AGB effects 
on top of these early models using the fuel consumption 
theorem. 

Besides the treatment of the TP-AGB, the BC03 and 
M05 models also differ for the set of stellar isochrones 
used to construct the synthetic populations, which may 
account for part of the difference in the results. 

In this paper we investigate further on the effect of stel- 
lar population recipes on the modeling of galaxy SEDs, 
and we do so by applying the two sets of models to infer 
the mass and age (or best-fit star-formation histories) to 
a specific sample of high-redshift galaxies. Note that the 
results obtained with the BC03 models are representative 
of what would be obtained with other widely used mod- 
els, like Pee ase llFioc fc Rocca-Volmerangd I1997D and 
Starburst99 i|Vazauez fc Leithereill2005D due to similar 
recipes for the TP-AGB phase and identical input stellar 
evolutionary tracks (M05). 

In Section 2 we recall the main differences between the 



BC03 and the M05 models, whereas Section 3 describes 
the main observational characteristics of the galaxies se- 
lected for this study. In Section 4 we describe the fit- 
ting procedure used to estimate the properties of these 
galaxies and the results are presented and discussed in 
Section 5. Finally, Section 6 is dedicated to a general 
discussion of the main results of this investigation. We 
adopt the current concordance cosmology, with Q\, Om 
and h (= i?o[kms-iMpc-VlOO]) equal to 0.7, 0.3 and 
0.7, respectively. The age of the best-fit model is re- 
quired to be lower than the age of the universe at the 
given spectroscopic redshift. 

2. MAIN DIFFERENCES BETWEEN THE TESTED STELLAR 
POPULATION MODELS 

The BC03 and the M05 stellar population models dif- 
fer by the following aspects: i) the stellar evolutionary 
models used to construct the isochrones; ii) the treat- 
ment of the TP-AGB phase; and iii) the procedure used 
for computing the integrated spectra. This latter point 
will not be discussed separately here (see M05). For full 
details on what below, we refer the reader to M05. 

Input stellar models. The B C03 models are based 
on the Padova stellar tracks fe.g.. iFagotto et an il994'l. 
whereas the M05 mode ls are mos tly based on the Fras- 
cati stellar tracks fe.g.. lCassisi et al. 1997') . The two set 
differ mostly for: a) the Padova tracks include a certain 
amount of convective overshooting on the Main Sequence 
(MS), whereas the Frascati tracks were constructed as- 
suming no overshooting; b) the temperature distribution 
of the red giant branch (RGB) phase, that is shifted to 
cooler temperatures in the Padova tracks for solar metal- 
licity and above, a result of the different calibration of the 
mixing-length theory used for the modeling of envelope 
convection. 

Convective-core overshooting has two effects on stellar 
models. First, it prolongs the MS lifetime for given stel- 
lar mass, while leaving the effective temperature almost 
unchanged. This results from the extended convective 
core ensuring more hydrogen fuel to the central regions 
where nuclear burning takes place. Thus, for given age 
of the MS turnoff, isochrones constructed from models 
with overshooting are hotter, therefore bluer, than those 
constructed from models without overshooting. Second, 
overshooting delays the appearance of stars with a degen- 
erate helium core, and the accompanying development of 
the RGB. For example, the RGB develops at an age of 
t ^ 0.5 Gyr in the Frascati m odels, and at an age of i ^ 1 
Gyr in the Padova models IjFerraro et alJ 12004 M05). 
The development of the RGB makes the SED redder, so 
in the age range ^ 0.5 — 1 Gyr models with overshooting 
are bluer, as they lack a well-developed RGB. 

The net effect is that theoretical SEDs based on over- 
shooting models are bluer therefore give systematically 
older ages compared to those based on models without 
it. Furthermore, by having a reduced contribution by the 
RGB, population models using stellar tracks with over- 
shooting have lower luminosity for given age and mass in 
stars, i.e., have higher M/L ratios, hence higher galaxy 
masses are derived for given luminosity. In summary, 
the net effect of overshooting is that for a given SED the 
BC03 models indicate older ages and higher masses com- 
pared to M05 models. Of course, this applies to popula- 
tion ages such that MS stars do have a convective core. 



i.e., for ages younger than a few Gyr. 

While it is hkely that overshooting makes convective 
cores somewhat more extended in mass, the size of this 
extension remains somewhat conjectural. On the one 
hand there exists no firm prediction from first princi- 
ples, like in most situations involving turbulence. For 
stellar populations between ^ 1 and ~ 10 Gyr, the cor- 
responding stellar masses at the MS turnoff range from 
~ 2 to ~ 1 Mq, with the mass of the convective core de- 
creasing with decreasing mass, and eventually vanishing 
for M r^ 1.2 Mq, the precise value depending on com- 
position. Thus, what one would need to know is the 
amount of convective overshooting (e.g., expressed as a 
fraction of the mass of the convective core) as a function 
of the stellar mass. This unknown function can only be 
empirically constrained by observations. In a recent at- 
tempt, based on fitting t he shape of the colo r -magn itude 
of some open clusters, iVandenBerg et all l|2006j) con- 
cluded that overshooting would amount to ~ 10% of 
a pressure scale height, of course vanishing as the con- 
vective core itself vanishes. An overshooting about this 
size i s also indicated by the asteroseismology of a /3 Gep 
star IjDupret et al.ll200^ . and by and large is consistent 
with the amount assumed in the Padova models. On the 
other hand, observations of LMC star clusters rule out 
the late RG B development implied by the Padova tracks 
IjFerraro et al. 2004) . To decide what is the size of over- 
shooting as a function of stellar mass goes well beyond 
the aims of this paper, and we shall limit ourself to ex- 
ploring the differences in the results when using different 
sets of models. 

The Temperature of the RGB. Whereas for ages 
less than ^ 1 Gyr the M05 code gives redder SEDs com- 
pared to the BC03 code, after the completion of the RGB 
phase transition this tendency is reversed and the M05 
SEDs are bluer (Figure 27 in M05). This is due to the 
warmer RGBs in the Frascati stellar tracks compared to 
the Padova ones (Figure 9 in M05), which in turn results 
from different calibrations of the mixing-length parame- 
ter used in modeling envelope convection. The ages de- 
rived with the M05 models are older than those obtained 
from the BCG3 models. 

The TP-AGB phase. The treatment of the TP-AGB 
phase is perhaps the main difference between the M05 
and the BC03 models. The contribution of the TP-AGB 
stars to the integrated light of a synthetic stellar popula- 
tion critically depends on what is adopted for the stellar 
mass loss taking place during this phase. The higher the 
assumed mass loss, the sooner the star loses its envelope, 
and the sooner the TP-AGB phase is terminated. Once 
more, there is no theory relating the mass loss rate to 
the basic stellar parameters, and once more one has to 
rely on empiri cal calibrations. BG0 3 adop t the mass loss 
calibration bv IVassiliadis fc WoodI ljl993|) . tuned to re- 
produce the maximum TP-AGB luminosity of a sample 
of Magellanic Cloud clusters from Frogcl ct al. (1990). 
In M05 instead, what is calibrated is the fractional con- 
tribution of the TP-AGB to the total bolometric light, 
and this calibration is made essentially on the same MC 
clusters using the Frogel et al. data. Sharing basically 
the same observational calibrating data, the two codes 
should give consistent results. As shown in M05 they 
do not, and the effective contribution of the TP-AGB is 
much larger in M05 compared to BC03. Full illustration 



of the differences is given in M05, where it is also shown 
that other population synthesis models, namely Pegase 
and Starburst99 behave very much like the BC03 ones. 
Again, rather than tracing the origin of the discrepancy, 
we focus here on the effects on the study of high redshift 
galaxies. 

Summing up the effects of both overshooting and TP- 
AGB, the M05 models are brighter and redder than the 
BC03 models for ages between ~ 0.2 and ~ 2 Gyr. Here 
the use of the M05 templates implies the derivation of 
lower ages and stellar masses. At older ages the M05- 
templates are instead bluer. 

3. GALAXY DATA AND SED FITTING 

For the present test we use the sample of seven galax- 
ies in the Hubble Ultra Deep Field (HUDF) which were 
singled out by the BzK criterion for z > 1.4 passively 
-evolving galaxies IjDaddi et alJl200^ . These galaxies 
show early-type morphologies on HUDF ima ges and have 
spectroscopic redshifts between 1.4 and 2.5 l|Daddi et alJ 
^05b, hereafter D05). Redshifts were determined from 
the Mg-I-Fe absorptions at rest-frame A -- 2600 - 2800 A 
producing a characteristic shape dubbed Mgjjv feature 
by DOS, on the HST-I-ACS grism spectra acquired in the 
framework of the Grism ACS Program for Extra-galactic 
Science project (GRAPES, Pirzkal et al. 2004). The 
presence of the Mgjjv feature imply that the UV/optical 
spectra of these galaxies are consistent with being domi- 
nated by A- or F-type stars. Therefore these objects are 
not dominated by OB type stars like typical star-forming 
galaxies and are predominantly passive. We choose these 
particular sample for this study because we might expect 
that the typical ages of stars in these galaxies are close 
to the critical 0.2-2 Gyr range where the different evo- 
lutionary synthesis models differ by the largest amount. 
This is therefore an ideal sample of high redshift galaxies 
to search for direct evidence of the effects from TP-AGB 
stars. 

Three of the seven galaxies were assigned a lower qual- 
ity ("B" class) redshift by D05, partly for the non perfect 
agreement with the photometric redshifts. This is now 
improved in most cases with the addition of the IRAC 
photometry to the SEDs and the M05 models (see Ta- 
ble 2 and Figure 9). 

Stellar population properties for these galaxies were 
derived in DOS by fitting their SEDs (ten bands from B 
to K) to BC03 templates spanning a wide range of pa- 
rameters. Reddening was restricted to be E(B—V) < 0.2 
and the Calzetti's law l)Calzetti et aU2 000') was adopted. 
In summary, the galaxies were found to be ~ 1 Gyr old, 
to be in passive evolution since at least ~ 0.4 Gyr, to 
have de Vauculeurs light profiles and large stellar masses 
( <; 10^^ -^0). They were therefore interpreted as progen- 
itors of massive ellipticals. The physical sizes of some of 
the objects {^ 0.7 kpc) were found to be significantly 
smaller than objects with similar masses in the local 
universe, which remained as an open issue in DOS. In 
the following we take advantage of the Spitzer IRAC 
(and MIPS) imaging obtained over the Ultra Deed Field 
(UDF) as a part of the Great Observatories Origins Deep 
Survey (GOODS) survey (Dickinson et al., in preparation 
and Chary et al., in preparation, respectively). All the 
seven galaxies were detected with high signal-to-noise in 
the four IRAC bands (see Figure QJ. IRAC photometry 
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Fig. 1.— The Spitzer IRAC and MIPS images of the seven HUDF 
galaxies. See Fig. 8 of DOS for HST ACS and NICMOS imaging. 



was obtained using 4" aperture magnitudes, corrected to 
total using aperture corrections for unresolved sources. 
The optical to near-IR SEDs of D05 were based on Sex- 
tractor MAG_AUTO apertures. We measured large 6" 
diameter magnitudes in the iiT-band to account for the 
flux lost by MAG_AUTO, to perform the match with the 
IRAC magnitudes. In order to account for the uncertain- 
ties in this matching process, and in the IRAC aperture 
corrections, we added a term of 0.1 mags in quadrature 
to the errors of the IRAC magnitudes. 

Figure |5] shows the observed I — K vs. K — A.hixm 
colour-colour diagram for our sample. Overplotted are 
some stellar population models from the M05, BC03 and 
Pegase libraries, with ages between 0.1 and 3.5 Gyr, in 
the observed frame at redshift 1.9. Simple Stellar Pop- 
ulations (SSPs) with solar metallicity are indicated by 
solid lines, while dashed lines display exponentially de- 
clining models with e-folding time of 0.3 Gyr (for the 
Pegase i nodels. only t he SSP is s hown). As in aiialogou s 
plots bv lLabbe etall l)2005j) and lPapovich et alJ l)2006j) . 
the very red K — 4.5^m colours displayed by the galaxies 
are not reached by the BC03 (or Pegase) templates, un- 
less one adds strong dust reddening to the models. The 
M05 models match instead the observed galaxy colours 
at the ages corresponding to the AGB phase-transition. 
The present database (see the Appendix for the photo- 
metric data) allows the sampling of the galaxy spectral 
energy distributions up to the rest-frame K-hand. The 
whole wavelength range from the rest-frame UV to K 
will be analysed in the next Section. 

3.1. SED fitting 
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Fig. 2. — Observed I — K vs. K — 4.5/im colour-colour diagram 
for the galaxies of our sample. Also shown are stellar population 
models from different libraries (M05, BC03 and Pegase) and for 
different star formation histories, placed at the illustrative redshift 
of 1.9, with ages ranging between 0.1 and 3.5 Gyr.. The vector 
indicates a reddening of A\j = 1 mag for a Calzetti et al. (2000) 
law ffrom lLabbe et alj|2005l) . 



We use an adapted version of the code Hyper-Z 
IjBolzonella et al.ll2000(l . kindly made available to us by 
M. Bolzonella, which performs SED fitting at fixed spec- 
troscopic redshift^. The fitting procedure is based on 
maximum-likelihood algorithms and the goodness of the 
fit is quantified via x^ statistics (see the Hyper-Z man- 
ual for details). The code computes the x^ for a certain 
number of templates, which differ for star formation (SF) 
histories, metallicities and ages, and finds the best-fitting 
template among them. It is important to note that the 
code does not interpolate on the template grids. This im- 
plies that the template set must be densely populated. 

We have constructed templates that cover a wide range 
of stellar population parameters, namely: i) metallicities 
from Zq/5 to 2 Zq; ii) ages from 10^ yr to 5 Gyr (the 
maximum possible age for the lowest redshift galaxy)^; 
iii) SF histories as: (a) instantaneous and chemically 
homogeneous bursts (i.e. Simple Stellar Populations, 
SSPs); (b) exponentially-decHning modes (SFR oc e"*/"^, 
or r-models, with r = 0.1, 0.3, 1 Gyr); (c) truncated mod- 
els, where SF is constant for a finite time interval ttrunc 
and zero thereafter (ttrunc = 0.1,0.3,1,2 Gyr); constant 
star formation ^. These composite templates (with pro- 
longed SF histories) were constructed using the BC03 
software package, for consistency with D05 and in order 
to be able to run the Hyper-Z code. As in D05, we deter- 
mined 95% confidence ranges for ages and masses follow- 
ing |Avni, (1976) . marginalizing over E(B-V), metallicity 
and SF history. This implies that the determination of 

^ For the fitting procedure we use photometric errors of 0.05 if 
the formal error is smaller than that, to account for systematics in 
photometry and color matching 

^ The spacing in age during the runs is fixed by the Hyper-Z 
code, see the manual for details. 

^ The models are available at www-astro. physics. ox. acuk/^ 
maraston 



the latter parameters is formally less constrained and the 
resulting values should be regarded as indicative. 

To allow a direct comparison with other works, DOS 
in particular, the templates refer to a straight Salpeter 
IMF down to 0.1 Mq. However, several pieces of ev- 
idence, from the redshift evolution of the fundamental 
plane to the M/L of local Es to the metal content of 
the intergalactic medium strongly favour a Salpeter-like 
IMF down to 1 Mq but flatter below as required to agree 
with the dynamically estimated M/L ratios of early- type 
galaxies t;Renz ini 2005). The actual shape of the IMF 
does not affect appreciably the SED since at any age most 
of the light is produced by stars in a narrow mass inter- 
val around the turnoff mass. However, it impacts on the 
amount of light emitted per unit mass turned into stars, 
thereby affecting the s tellar M * IL ra tio. For example, 
when using the IMF of lKroupal l|2001j) the M/L ratio is 
a factor ~ 1.5 smaller than in the straight-Salpeter case 
(M05). 

For the objects classified of 'Class A redshift' by D05 
the fitting procedure included also the strength of the 
Mguv absorption feature, as defined in D05 and mea- 
sured on the GRAPES spectra. The difference between 
observed and model Mgyy strength was included in the 
calculation of the x? when searching for the best-fit 
model, with the same weight as the photometric data. 

Finally, the entire SED fitting was repeated with 
BC03-based templates spanning the same stellar popu- 
lation parameters as those used for the M05 models. 

The output of the procedure includes: the age t, i.e., 
the time since the start of SF; the metallicity [Z/H]; the 
star-formation history SFR; the reddening E{B—V)\ and 
the stellar mass M* . 

The stellar mass M* is evaluated by comparing the ob- 
served and the synthetic SED, which provides the mass 
that went into stars by the age of the galaxy. This ap- 
proach was retain in order to allow a comparison with 
D05. However, such mass overestimates the true stel- 
lar mass, as stars die leaving remnants whose mass is 
smaller than the initial one. Adopting the prescriptions 
by Renzini & Ciotti (1993) for the initial mass-remnant 
mass relation, M98 evaluated that the stellar mass of a 
SSP is reduced by ^ 30% 15 Gyr after formation, most 
of which happens within the first few Gyrs. For given 
mass turned into stars, such reduction is smaller in stel- 
lar populations with extended star formation histories, 
due to a relatively higher fraction of living stars. The ac- 
tual M*{t) for some illustrative cases of star-formation 
histories is given in Figure |21 Tables 1 and 2 give the 
stellar mass prior of such reduction, thus enabling a di- 
rect comparison with other results. Table 2 also reports 
the mass decrement (in percent) that corresponds to the 
given SFH and age^. This mass decrement should be 
applied to the masses listed in column 12. 

4. RESULTS 

4.1. Assuming no reddening, E{B — V) — Q 

We first discuss the solutions that are obtained in the 
assumption of zero reddening. This allows us to single 
out the effect of stellar population modeling, thus high- 
lighting some interesting differences between the two sets 

** Tabular values of M* for several star formation histories can 
be found at www-astro. physics. ox. acuk/f^ maraston. 
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Fig. 3. — The actual stellar mass M* , normalised to an initial 
value of 1 Mq , of stellar populations with different star formation 
histories. 



of models. The corresponding best-fit solutions (lowest 
X^) for using the BC03 (blue fines) and the M05 (red 
line) models are shown in Figure 0] The observed SEDs 
are shown as solid symbols with errorbars, and informa- 
tion on the individual galaxies and best fit SFHs is given 
in Table 1. Overall, the M05 models allow for better 
fits (lower Xr) for all galaxies, in the rest-frame near- 
IR as well as in the blue (e.g., the objects at z=1.39, 
1.55 and 2.47). Given the ages required to fit the ob- 
served SEDs, it is likely that the different recipes for 
the TP-AGB phase are primarily responsible for such ef- 
fect. For the objects at z = 1.91 and z = 2.47 none of 
the BC03-based templates was able to match the near-IR 
fluxes. In the case of the z = 2.47 galaxy not even an age 
as old as 2.3 Gyr is sufficient to obtain enough near-IR 
flux. Note also that the SSP solution gives a fairly good 
flt with M05 models (which are TP-AGB dominated at 
this age), whereas an extended SF history (constant un- 
til itrunc = 2 Gyr) was required when using the BC03 
models. In this case the M05 single-burst models supply 
both the high optical and the high near-IR fluxes, which 
come from the warm turnoff and cool TP-AGB stars, re- 
spectively. The stellar masses M* are typically smaller 
for the M05 models, being on average ~ 60% of those 
derived with the BC03 models. In part this is a conse- 
quence of the lower ages, with the average age of M05 
models being also ~ 60% of that derived for the BC03 
models. The differences are not as high as they would 
be in case of single bursts, as the more extended SFHs 
required by the BC03 models have lower M* /L ratios. 

The most relevant difference in the results from the two 
sets of models is about the implied formation redshifts 
(see Table 1). Using BC03 models the stellar popula- 
tions form at higher redshifts (sometime at much higher 
redshifts, see e.g., the object at z = 2.47) and, since the 
galaxies are also more massive compared to the case when 
using M05 models, an appreciably faster mass growth of 
galaxies is implied. 
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Fig. 4. — The observed spectral energy distributions (filled symbols with errorbars) and the best fit solutions obtained with the M05 
and the BC03 templates (red and blue lines, respectively). The parameters of the fits are indicated on the figure and given in Table 1. No 
reddening was adopted (i.e. E(B-V)=0). 



4.2. E{B ~ V) as a free parameter 

The best-fit solutions when also reddening is treated 
as a free parameter are displayed in Figure [S] and the 
parameters of the best-fits are given in Table 2. The 
reddening was allowed to vary among the options offered 
in the Hyper- Z package, namely, Milky Way, SMC, LMC 
and the so-called Calzetti's law for local starbursts (refer- 
ences in Table 2). All these prescriptions refer to a dust- 
screen configuration. For this 'reddening' case, in order 



to exploit all the available information and to control the 
degeneracies introduced by the larger number of free pa- 
rameters, we have additionally included the strength of 
the Mg(jv feature in the computation of the x^ ■ This 
was done only for galaxies with Class A redshift, namely 
for objects 4950, 1025, 3650, 3574, because for the Class 
B objects the Mguv is a less robust measurement (D05). 
By allowing the best-fit procedure to treat E{B — V) as 
an additional free parameter, and allowing it to choose 
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among four different reddening laws, it is no surprise that 
solutions with appreciably better x^ £^re found. However, 
this does not ensure that such solutions are more likely 
than those assuming E{B — V) — 0. It is also no surprise 
that some best-fits require the Calzetti's law, whereas 
others prefer the SMC law. In just one case the LMC law 
was preferred, while in no cases the MW law was chosen. 
It is worth noting that the best-fit for the z = 2.47- 
object exhibiting the strongest TP-AGB contribution is 
reddening-independent for the M05 models. 



In general, the inclusion of reddening reduces the ages 
of the best-fits obtained with the BC03 models, mak- 
ing them closer to those obtained with the M05 models. 
However, reddening does not help to improve substan- 
tially the fits to the IRAC bands, like in the case of ob- 
jects 1446 and 3650. TP-AGB effects cannot be traded 
with dust or metallicity effects, which helps reducing the 
degeneracies in the age determination. Detailed com- 
ments on the fits for individual galaxies are given in the 
Appendix. 
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Fig. 6.— The GRAPES HST+ACS spectrum of galaxy 3650 
at 2 ^ 1.91 (from IDaddi et alj 12005bl) 'l. with superimposed the 
best-fits obtained with the BC03 and the M05 models (blue and 
red lines, respectively, smoothed at the same spectral resolution of 
the observed spectrum). 



Figure El compares the GRAPES HST+ACS spectrum 
of object 3650 at z ~ 1.91, with the best-fit solutions 
corresponding to the BC03 and M05 models (see Ta- 
ble 2). We use this plot for illustrating the effect of the 
age/metallicity degeneracy that affects the optical, but 
is alleviated in the near-IR due to the distinctive spec- 
tral features of TP-AGB stars. The best-fit solutions 
recover similarly well the strength of the Mgyy line and 
the spectral shape around it, in spite of having different 
population parameters, namely young and metal-rich for 
the M05 models and older and metal-poor for the BC03 
models. However, the BC03-based solution gives an ap- 
preciably poorer fit to the near-IR SED (cf. Figure O, 
showing how the inclusion of the near-IR fluxes can help 
breaking the age/metallicity degeneracy. 

However in most cases spectral absorptions and pho- 
tometric data on a wide spectral range are not available 
simultaneously. How realistic are then the best-fits ob- 
tained from the sole photometry? To get insight on this 
issue, we have repeated the fitting procedure without in- 
cluding the Mguv line when searching the best-fit model. 
We have then checked how well these 'photometric' best- 
fit models recover the observed Mguv ■ This is shown 
in Figure in which the Mgyy strengths predicted by 
the best-fit 'photometric' models are compared to the ob- 
served ones. Although based on a very small number of 
objects, the predicted values from M05 models are some- 
what more accurate than those from the BC03 models 
(variances are 0.04 and 0.08, respectively). 

As a final exercise, we have checked the effect on photo- 
metric redshifts derived using one or the other set of pop- 
ulation models, and of the inclusion of the IRAC fluxes. 
The results are given in Table 2 and shown in Figure |S1 
The photometric redshifts obtained with the M05 tem- 
plates are in better agreement with spectroscopic red- 
shifts than those obtained with the BC03 models. The 
photometric redshifts obtained with the BC03 models 
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Fig. 8. — Comparison between the spectroscopic redshifts and 
the photometric redshifts obtained with the BC03 and M05 models. 
The IRAC magnitudes have been included in the derivation of the 
photometric redshifts. 



and the IRAC data appear to be more discrepant than 
those in which only magnitudes up to K were used (cf. 
D05). Also this effect likely originates from the different 
TP-AGB recipes. 

4.3. Stellar Masses 

The different star-formation histories derived from the 
different models have an impact on the derived stellar 
masses, which is summarised in FigureEl The upper pan- 
els display the effects of the models (i.e., M05 vs BC03) 
on the derived stellar mass of galaxies, when using only 
the optical and near-IR bands (left-hand panel) and when 
also the IRAC 3.5-8 /zm data are used (right-hand panel). 
The masses predicted by the M05 templates are lower, 
which is a consequence of the lower ages. The effect is 
appreciably more pronounced when the IRAC bands are 
included, as the TP-AGB phase (the most discrepant in- 
gredient between the two models) has its largest contri- 
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bution in the rest-frame near-IR. However, lower masses 
can also be derived when the observed-frame K is the 
reddest band used in the best fit. In fact, for the objects 
at z = 1.39.1.55, 1.73 the K band samples the rest-frame 
I, which is aheady affected by the TP-AGB (M05, Fig- 
ure 18). In addition, the later onset of the RGB phase in 
the overshooting tracks adopted by BC03 helps to make 
the BC03 templates fainter at ages around 0.5 Gyr. The 
corresponding M/L ratio is higher and a higher M* is 
derived (see Section 2). It is embarrassing to realize to 
which extent the recipes for the convective overshooting 
impact on the stellar masses that are derived for high-z 
galaxies, and therefore on the implied galaxy formation 
scenarios. 

The lower panels in FigureElshow the effects (for given 
stellar population models) of the wavelength range in- 
cluded in the fitting procedure: only up to the K band 



(on the X axis) and with the IRAC bands (on the y axis) . 
The results based on BC03 are not appreciably affected 
by the inclusion of the IRAC bands (as recently pointed 
out bv lShaplev ct al. 2005). Instead, for the M05-based 
solutions, three galaxies arc found to have appreciably 
lower masses when the IRAC fluxes are considered, again 
the likely effect of the inclusion of the TP-AGB contri- 
bution. 

Finally, Figurc^jcompares the masses obtained in this 
paper (including the IRAC data) with the c alibration for 
z ^ 2 BzK-selected galaxies proposed by iDaddi et alJ 
l)2004l) . This calibration was obtained using the observed 
-ftT-band magnitude and z — K colour only, along with the 
BC03 templates. The figure shows that the calibration 
is well consistent with the BC03-based masses derived in 
this paper, which is no surprise. The M05 models suggest 
instead a calibration lower by 0.2 dex, on average. 
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4.4. MIPS detection and ongoing SFR 

Only two galaxies are significantly detected with 
Spitzer+MIPS at 24^m (objects 1025 and 4950, Fig.[TJ. 
The MIPS 24^m fluxes of about 25-30/iJy convert 
to star-formation rates of about 20-30 Moyr"^ at 
the redshifts of the sources, using the models from 
ICharv fc Elbaz (2001) . The two MIPS detections are in 
fact the galaxies with the bluer UV spectra (see D05) 
and more extended star formation histories (see Ta- 
ble 2). The MIPS derived SFRs are comparable, given 
the uncertainties, with the residual ongoing SFRs (of 
~ 10 - 2OM0/yr) that are obtained through the SED 
modeling of these 2 galaxies (see also D05). For the rest 
of the galaxies the 3 a limits on the 2Afj,m flux imply SFR 
below 10 MQyi~^. We recall that typical z — 2 galaxies 
with similar K-band fluxes to the sample studied here 
have much stronger 24/z?n emission, of order of 120/iJy 
l|Daddi et al. 2005ai). The low MIPS fluxes from these 
sources underline once again, in agreement with their 
spectral features, that these objects are predominantly 
passive. 

4.5. The role of AGNs 

In obtaining the best-fltting models we have assumed 
that the observed SEDs are the results of just stars and 
dust. The possible presence of an AGN would obviously 
affect stellar population modeling in general. However, 
differences between the physical properties of galaxies 
inferred with the M05 or BC03 models remain, as these 
models have a different ratio of optical to near-IR flux. 
Therefore, even if at play, AGNs would not alter the sense 
of this investigation. 

More specifically, AGNs probably play a minor role on 
our analysis. First, no clear AGN features are observed 
in the spectra of these objects. Also, the shape of the 
SED over IRAC is clearly stellar with a secure identifica- 
tion of the peak emission at 1.6 /im. Two of the galaxies 
(1025 and 1446) have AT-ray detections, but as discussed 



in D05 these are extremely hard and absorbed. Cou- 
pled to the evidence of the faintness of the galaxies at 
24 fim, we conclude that the AGN should not influence 
appreciably the photometry from the UV to the near-IR 
rest-frame. For example, if we adopt the SED shapes of 
NGC1068 or MRK231 and normahse them to the 24^m 
photometry of the sources, the contribution of the AGN 
to the optical to near-IR light of the galaxies would be 
entirely negligible. Similarly, from the nar row JtT-band 
to X- ray flux correlation of AGN shown by iBrusa et alJ 
()2005|) . we conclude that possible AGN in our sources 
are expected to provide a negligible contribution to the 
ivT-band light of our galaxies. 

Finally, the overall good agreement between predicted 
and observed Mgyy line-strengths (Figure 7) and the 
good match between spectroscopic and photometric red- 
shifts (Figure 8) suggest that a significant AGN contri- 
bution is highly unlikely. 

5. DISCUSSION 

The aim of this paper is to explore how much the star 
formation histories of high-z galaxies, that are derived by 
fitting observed to synthetic spectral energy distributions 
depend on the adopted stellar population models. In par- 
ticular, we focused on the TP-AGB evolutionary phase, 
because TP-AGB stars are the dominant bolometric and 
near-IR contributors in stellar populations with ages ~ 
1 Gyr (M98; M05). Therefore, such phase must be of 
primary importance at hig h-z, when the universe itself 
was j ust a few Gyrs old (Maraston ' 2004 l2005t IRenziii] 
11992^ . and especially in the interpretation of Spitzer data 
for high-z objects. The recipes for the TP-AGB phase 
are found to differ significantly among the various mod- 
els, which has motivated us to check the robustness of 
the results obtained for high-z galaxies. 

To this aim we have analysed a sample of mostly pas- 
sively evolving high-z (1.4 ^z;^ 2.5) galaxies with opti- 
cal and Spitzer/IRAC photometry, modeling the whole 
spectral energy distribution from the rest-frame B to K . 
The analysis was performed with two libraries of stellar 
models: the M05 models in which the TP-AGB phase 
is calibrated on Magellanic Cloud globular clusters, and 
the BC03 models in which, like in most other models in 
the literature (e.g. Pegase, Starburst99) the energy con- 
tribution from the TP-AGB phase is significantly lower 
than in the Maraston models. Furthermore, the two 
models adopt stellar evolutionary tracks, in which the 
onset of the Red Giant Branch is predicted at different 
ages, which also impact on the time evolution of the syn- 
thetic spectral energy distribution. The differences be- 
tween stellar population models are widely discussed in 
M05 and recalled in Section 2. 

Besides the 14 photometric bands spanning from rest- 
frame B to K , the procedure introduced in this paper 
also includes the strength of a near-UV absorption fea- 
ture (the Mgyy at 2800 A, see D05) for searching the 
best-fit galaxy models. In this way the effect of the 
age/metallicity degeneracy is substantially reduced. We 
have also checked that even without using the spectro- 
scopic information, but only with the photometric SED, 
the best-fits obtained with the M05 models recover well 
the strength of the Mguv feature, hence provide consis- 
tent models for distant galaxies. 

We then show to which extent the derived star- 
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formation histories of the sample galaxies are model de- 
pendent. The M05 models typically match the rest-frame 
near-IR and optical very well, without invoking strong 
reddening and without violating the constraint set by 
the age of the universe. The BC03-based solutions ex- 
hibit in some cases near-IR fluxes that are significantly 
lower than the data, which requires the addition of old 
stellar populations or dust reddening or a combination 
of the two. This is reminiscent of recent results on high- 
z galaxies in which Spitz er/IRAC data a re interpreted 
with the BC models (e.g..lYan et al. l l20M) or the Pegase 
models fe.g.. lVillar-Martin et al. 200(|). The young solu- 
tions required to match the blue side of the SED could 
not match the rest-frame near-IR, which led the authors 
to assume an underlying old population to be responsible 
for the near-IR fluxes. While a two-component model is 
not implausible in itself, the proportions by mass of the 
young and old component (1 vs. 99 %) as well as the 
age of the old component, comparable to the age of the 
Universe, make the star formation history contrived. 

The lower masses derived from the M05 models may 
help understanding the origin of an apparent discrepancy 
noted by DOS and concerning the relationship between 
the physical sizes and the masses of the galaxies in the 
present sample. The galaxies were found to be too small 
for the given masses in comparison to local objects. The 
lower masses derived here already help in alleviating the 
discrepancy. In particular, the z ~ 2A7 object, with 
a stellar mass around 10^° Mq, is more likely the pre- 
cursor of a bulge rather than a massive elliptical. Its 
rather small size {r^ ~ 0.8 kpc, D05) supports this con- 
clusion, as it can be inferred from comparisons with local 
samples IjBender et al.lll99^ . Also its formation redshift 
and metallicity compares well with similar estimates of 
bulges in the local universe IjSaraiedini fc Jablo nka'20051 
IXhomas & Davics 2006). According to the BC03 models, 
three galaxies have rather low metallicities (Zq/5) and 
at the same time high masses (> 10^^ Mq). Such objects 
do not match local relationships, even taking metallicity 
gradients into account. Only one such object is found 
with the MG5 models. 

Finally, we find that the stellar population templates 
impact on photometric redshift determination. The use 



of IRAC bands and M05 templates give photometric red- 
shifts that are in excellent agreement with the spectro- 
scopic ones. Again, we ascribe this effect to the prescrip- 
tions for the TP-AGB phase. 

In this work we have focused on high- 2 galaxies, where 
the stellar populations are young therefore TP-AGB 
stars must exist. Our spectral fitting represents the first 
direct evidence for such stars in the early universe. Go- 
ing to even higher redshifts might reveal further inter- 
esting aspects of galaxy evolution, a scope we are pur- 
suing in the future. However, TP-AGB stars may mat- 
ter at any redshift depending on the specific star for- 
mation history. For example, Ivan der Wei et all |)200,'tI ^ 
measured the evolution of the rest-frame i^-band Fun- 
damental Plane from z ~ 1 to the present, by using 
IRAC imaging of a sample of early- type galaxies. They 
find that the K-hand luminosities of the low-mass ob- 
jects are too large (at given B luminosity), or the B — K 
colours t oo red, with r e spect to the prediction of the 
BC03 or lVazdekisei~ai] ((l99l models. The M05 mod- 
els were found to provide a better match the data. This 
effect, which is correctly attributed to the treatment of 
the TP-AGB phase, implies that those low-mass early- 
type galaxies experienced a burst of star formation at 
redshift close to one. 

In conclusion, the recipe for the stellar population 
modeling plays a central part in the derivation of the 
star formation histories of high-z galaxies, which should 
be kept in mind when the derived parameters are used 
to constrain galaxy formation. 
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Table 1. Galaxy SED fitting results: no reddening 



ID'' 



mod 



^for 



Gyr 



[Z/H] 



SFH 



M0/IO" 



M.t 



8238 1.39 



4950 1.55 



1025 1.73 



3523 1.76 



3650 1.91 



3574 1.98 



1446 2.47 



M05 
BC03 
M05 


1.0 
2.3 
1.7 


1.8 
2.9 
2.7 


2Z0 
2.5 Z0 
Z0/2 


g-t/O.l Gyr 
g-t/0.3 Gyr 
g-t/0.3 Gyr 


1.5 
2.6 
1.5 


U-40_o.o6 

2.20«i 


BC03 

M05 


1.7 
1.7 


2.7 
3.2 


2.5 Zq 

z© 


g-t/0.3 Gyr 
g-t/0.3 Gyr 


4.9 
2.2 


i.oo«-i 


BC03 
M05 
BC03 


1.7 
1.4 
2.3 


3.2 
2.9 
4.6 


2.5 Z0 

Z0 
2.5 Z0 


g-t/0.3 Gyr 

itrunc = 1.0 Gyr 
itrunc = 2.0 Gyr 


3.8 
2.2 
3.1 


1.25+J-f 
0.60|:|6 


M05 


0.5 


2.2 


2 Z0 


itrunc = 0.1 Gyr 


2.5 


'J-oO-o.oi 


BC03 


2.3 


5.5 


2.5 Z0 


ttrunc = 2.0 Gyr 


3.5 


^■"-0.5 


M05 


2.3 


6 


Z0/5 


itrunc = 1.0 Gyr 


2.6 


0.56t«:i^ 


BC03 


2.6 


8.2 


0.4 Z0 


ttrunc = 2.0 Gyr 


4.0 


0.53+Jo:\l 


M05 
BC03 


0.4 
2.3 


2.8 
15 


2Z0 
Z0 


SSP 
ttrunc = 2.0 Gyr 


6.2 
10.4 


0.60in^ 



Note. — The variable t is the temporal distance from the beginning of star-formation. The formation redshift 2:form I'cfcrs to t, thereby indicating, 
in case of extended star formation histories, the epoch at which the galaxy started to form. The last IRAC filter was excluded for objects up to 
z=1.76, for which this filter samples rest-frame wavelengths > 2.5 microns, that is the limit where the empirical TP-AGB star spectra are defined 
in the M05 models. The same filter configuration was used when the BC03 models were employed. We use the following conversion between AB 
and Vega magnitudes in the B- and iC-band: Mgi-vcga) — ^b(Ab) ~ ( — 0.089), A/j^fVcga) — ^^k(AB) ~ (1.871) . 

^Galaxy ID as in Daddi et al. 2005 

^The mass is obtained from normalisation of the SED. Ranges refer to the 95% confidence level. Note that the total magnitudes are 
defined from the values obtained from Sersic profile fitting given in D05, Table 1. 



Table 2. Galaxy SED fitting results 



ID'' 


^spcc 


■S^phot 


mod 


t 
Gyr 


^form 


[Z/H] 




SFH 


E(B-V) 


M0/IO" 


Af.<= 
M0/IOII 


M0/IOI1 


M,g 


%M^ 


8238 


1.39 


1.39 


M05 


0.7 


1.7 


2Z0 


ttrunc = 0.3 Gyr 


0.15*^ 


1.1 


0.25-0.47 


Q 45 + 0.22 


1-2.4 


21 






1.30 


BC03 


2.3 


2.9 


2.5 Z0 


itrunc = 2.0 Gyr 


0.15<^ 


2. 




1 n+0.45 
^•'J-0.40 






4950 


1.55 


1.55 


M05 


3.5 


7.5 


Z©/2 




g-t/lGyr 


0.15*^ 


0.6 


2.8-5.6 


5.00+»« 


2.9-7.3 


26 






1.67 


BC03 


3.5 


7.5 


Z0/5 




g-t/lGyr 


0.25^^ 


0.94 




8.10l»1 






1025 


1.73 


1.72 
1.77 


M05 
BC03 


1.4 
1.4 


3.1 
3.1 


Z0 
Z0/2.5 




g-t/0.3Gyr 
g-t/0.3Gyr 


0.07*^ 
0.15^^ 


1.1 
0.64 


0.3-0.4 


0.62lg:?8 

J^-J^-0.26 


1-1.8 


20 


3523 


1.76 


1.68 
1.63 


M05 
BC03 


0.3 
0.5 


1.9 
2.1 


Z0/5 
Z0 


itrunc = 0.1 Gyr 

g-t/O.lGyr 


0.4*= 
0.4"^ 


1.2 
1.2 


0.4-0.6 


om+°f 

0.92+J,it 


1-1.5 


18 


3650 


1.91 


1.88 


M05 


0.3 


2.1 


2 Z0 




SSP 


0.2*= 


1.3 


0.33-0.85 


0.83l°:« 


1.3-2 


19 






1.81 


BC03 


1.0 


2.7 


Z0/5 




SSP 


0.2'= 


1.9 




n O+0.14 

^•'J-o.gs 






3574 


1.98 


1.94 


M05 


0.2 


2.1 


2Z0 




SSP 


0.3'' 


2.5 


0.2-0.4 


„ 42+0.28 


0.5-0.9 


18 






1.94 


BC03 


1.0 


2.9 


Z0/5 


itrunc = 0.3 Gyr 


0.2== 


2. 




n 7+0.3 
^■'-0.2 






1446 


2.47 


2.74 
2.86 


M05 
BC03 


0.4 
1.0 


2.8 
3.8 


2 Z0 

Z0 




SSP 

g-t/O.lGyr 


0.00 
0.35= 


3.r 

4.7' 


0.27-0.2 


'^•^'-0.04 

1 n+o.4s 


0.7-1.1 


20 



Note. — The best-fit parameters and Xr refer to the model obtained using the spectroscopic redshift ^spcc- The values 2:phot arc given for 
comparison. See Tabic 1 for other notes. 

'^Galaxy ID as in Daddi et al. 2005 

'^SMC reddening law as in lPrevot eFall 119841) and lBouchet et al.l IT981) 

"^Reddening law as in Calzetti et al. (2000) 

•^LMC reddening law as in lFitzpatrlc^ 119861) 

<=Through Af*/LBv„g, and M*/L^^^^^ 

^See note (b) to Table 1 

sProm Daddi et al. 2005 

'^This column gives the percentage of mass decrement due to stellar mass loss. It should be applied to values in Column 11 (labelled f) 

'For this object the x^ were computed without the photometric data in the Nicmos i^llOW filter. This point is off but has a small 
errorbar, which drives the x^ to the misleading values of 6.2 and 5.8, in case of the M05 and BC03, respectively. 



APPENDIX 
Here we comnient in more details on the best-fits obtained for individual objects for the reddening case (cfr. Table 2). 
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8238, z=1.39. The Mgyy line was not included in the x^ The object is ~ 0.7 Gyr old, metal-rich and in passive 
evolution since 0.4 Gyr. The reddening follows an SMC-like law. The stellar mass is O.45TO^^M0 with small dispersion. 
For the BC03 models, the best solution has a worse x^ ; a significantly older age and a significantly longer formation 
timescale. The stellar mass is a factor 2.2 larger and consistent with what obtained in D05 from fitting the sole optical 
bands. 

It is interesting to note that the same, smaller stellar mass would have been derived by fitting the M05 models up the 
sole observed- frame K (see Figure Ej). This is due to two effects. First, at redshift ~ 1.4 the observed K samples into 
the rest-frame /, which is already affected by the TP-AGB prescriptions (M05). Second, as pointed out in Section 2, 
another important difference between the two sets of models exist, namely the energetics of the different input stellar 
evolutionary tracks. The M05 models are based on tracks in which the RGB phase starts developing at younger ages 
with respect to the Padova tracks used in BC03. The delay in the latter is due to the overshooting (see discussion in 
M05). The result is that the M05 SSP models are brighter around 1 Gyr (cf. M05, Fig. 7), which implies a lower stellar 
mass, an effect anticipated in M05. It should be noted that the earlier on set of the RG13 pha se compares very well 
with the observed RGB contribution in Magellanic Cloud globular clusters IjFerraro et alJl20M) . while the predictions 
based on the Padova tracks display a severe discrepancy with the data. 

4950, z=1.55. The Mgyy line was included in the x^ The best-fit solution is 3.5 Gyr old, is forming stars according to 
a tau model with e-folding time of 1 Gyr, the metallicity is half-solar and E{B — V) = 0.15 (SMC law). The extended 
star formation history is consistent with the irregularities and blobs in the _B-band image and with its spiral-like 
morphology (D05). This galaxy has a huge De Vauculeurs bulge and turned out to be the most massive object of the 
sample. The observed Mgyy is remarkably well reproduced by this model. A'L^ is large and very well constrained. It is 
consistent with both the value of D05 and that derived here with the BC03 models. The best solution obtained with 
the BC03 models has similar overall parameters, but the reddening is larger and the x^ is worse. 
1025, z=1.73. The Mguv line was included in the y^ The best-fit is ~ 1.4 Gyr old, still forming stars at a very 
modest rate (t model, with r = 0.3 Gyr). The stellar mass is 0.6 IO^^Mq. The solution obtained with the BC03 
model has a lower metallicity and a larger mass. 

3523, z=1.76. The Mgyy line was not included in the x^ The Mgyy for this galaxy is very smaU (0.094 ± 0.09 A), 
which implies either very young ages or very low metallicities. The solution allowing for the best comparison to the 
observed line turned out to be young and rather metal-poor. Its Mg^y is 1.07. The galaxy is passively evolving since 
0.2 Gyr. The solution obtained with the BC03 models allows a worse match to the Mgyy line (1.14 A) because of the 
higher metallicity. A higher metallicity is found because it provides more flux to the near-IR. 

3650, z=1.91. The Mguy line was included in the x^ The galaxy is young and metal-rich, and in passive evolution 
for the M05 models. It has formed stars at a remarkable rate as M* is 0.83 • 10^^ Mq. The spectrum displays a nice 
near-IR excess that is very well fitted. The match to Mguy is remarkable good (1.25 A versus the observed 1.25 A). 
The solution obtained with the BC03 models is older and significantly more metal-poor, though the stellar mass of 
the galaxy if higher (2 • 10^^ Mq). This mass is too large for a metallicity Z^/b in comparison with what found in the 
local universe. This solution matches worse to the observed SED. 

3574, z=1.98. The Mgyy line was included in the x^ The best-fit is young, metal-rich, and rather highly reddened. 
The model is passively evolving. The best-solution obtained with the BC03 models is older and metal-poor. 
1446, z=2.47. The Mgyy line was not included in the x^ This is perhaps the strongest case of AGB-dominated object. 
The fit requires a super-solar metallicity. However, the amount of TP-AGB in the M05 models at this metallicity was 
even slightly too low around the J band with respect to the observations. Since the TP-AGB fuel and spectral- types at 
super-solar metallicities are not well constrained given the absence of local GCs as calibrators (see discussion in M05), 
we could accept to modify the model. An excellent fit is obtained when the fuel is augmented by 30% and the 0-rich 
stars have a slightly warmer temperature. Both changes are certainly inside the uncertainties of the TP-AGB phase 
at high- metallicity Noteworthy, this fit does not require any reddening. The stellar mass is less than 1/3 the lower 
bound of the D05 solution. The BC03-based solution is older and highly reddened, yet the match to the observed SED 
is not satisfactory, in the near-IR as well as in the optical. We have verified that older ages (2.6 Gyr) would provide a 
better fit, but these ages are above the limit set by the age of the Universe. 

Note that the visual quality of the fit did not correspond to the x^ values when the FllO band was included in 
the x^ This is due to this photometric data having a very small error, which drives the x^ While the best-fit model 
was found using all bands for consistency with the other galaxies, the x^ values given in Table 2 were computed by 
excluding the FllO band. 

The data in FllO and J bands are off the models and are responsible for the significantly higher values of x^ that 
are obtained in case of this galaxy with respect to the other objects. We have noted that significantly better values of 
X^ (~ 1.5) are obtained when the redshift is higher than the spectroscopic value and close to the redshift estimated 
photometrically with the M05 models (2.74, cf. Figure IHJ. In this case, the best-fit model, that has similar overall 
parameters as that presented in Figure [S] did not require any modification of the fuel consumption in the TP-AGB 
phase. 
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Table 3. Photometric data of HUDF galaxies 



ID 


F435 


V 


F606 


R 


F775 


F895 


FllO 


8238 


29.15 ±0.62 


28.59 ±1.10 


27.42 ±0.09 


26.86 ±0.2 


25.60 ±0.02 


24.66 ±0.01 


23.67 ±0.01 


4650 


25.30 ± 0.04 


24.83 ± 0.04 


24.78 ±0.02 


24.47 ± 0.03 


23.89 ± 0.01 


23.18 ±0.01 


99 


1025 


26.35 ± 0.05 


25.66 ± 0.20 


25.84 ±0.02 


25.56 ±0.09 


25.17 ±0.01 


24.52 ±0.01 


99 


3523 


28.40 ±0.39 


28.33 ± 0.88 


26.77 ±0.06 


26.52 ±0.14 


25.73 ± 0.03 


24.99 ±0.02 


23.60 ±0.01 


3650 


28.51 ± 0.31 


27.00 ± 0.27 


26.42 ± 0.03 


26.10 ±0.10 


25.48 ±0.01 


24.57 ±0.01 


23.36 ±0.01 


3574 


28.81 ± 0.34 


28.49 ±1.04 


27.10 ±0.05 


27.31 ±0.29 


26.38 ± 0.03 


25.57 ±0.02 


24.41 ±0.01 


1446 


29.92 ±1.10 


28.20 ±0.78 


27.52 ±0.08 


26.63 ±0.12 


26.44 ±0.03 


26.12 ±0.04 


25.43 ± 0.03 



Table 3. Continued 



J _F160iy K m(3.6/im) m(4.5^m) m(5.8Atm) m{8.0^m) MIPS(24^m) 

23.16 ±0.04 22.76 ±0.0 22.08 ± 0.02 21.29 ±0.13 21.38 ±0.13 21.71 ±0.13 22.06 ±0.15 4.9 ± 3 

21.57 ±0.02 99 20.54 ±0.11 19.89 ±0.11 19.84 ±0.11 20.11 ±0.11 20.58 ±0.11 26.7 ±4.8 

22.74 ±0.07 99 21.80 ±0.05 21.28 ±0.12 21.20 ±0.12 21.27 ±0.12 21.62 ±0.13 31.4 ±3.5 

23.19 ±0.04 22.68 ±0.01 22.33 ± 0.03 21.60 ±0.14 21.56 ±0.14 21.69 ± 0.15 21.88 ±0.15 6.5 ± 5.8 

22.66 ±0.02 22.28 ±0.00 21.75 ±0.02 21.05 ±0.12 21.01 ±0.12 21.10 ±0.12 21.65 ±0.12 5.8 ± 4.3 

23.63 ±0.05 23.22 ±0.01 22.67 ± 0.03 22.24 ±0.15 22.19 ±0.15 22.12 ± 0.15 22.73 ±0.18 15 ± 6 

25.25 ±0.22 23.38 ±0.01 22.86 ± 0.04 22.18 ±0.19 22.08 ±0.19 22.01 ± 0.19 22.03 ±0.19 3 ± 4.6 

-"^■The magnitudes arc in AB system, which arc related to flux density fu (in erg s^-*^ cin^^ Hz^-*^) by m = —2.5 X lg{fi/) — 48.60. Fluxes 
at 24 /I are in /ijj/ 
^The reported photometric errors of the IRAC bands reflect the random errors only. Typical systematic errors in these bands arc at < 0.1 



mag level. 



